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O i We introduce four types of SU(2M + 1) spin chains which can be regarded as the 

BCn versions of the celebrated Haldane-Shastry chain. These chains depend on 
two free parameters and, unhke the original Haldane-Shastry chain, their sites need 
0^ I not be equally spaced. We prove that all four chains are solvable by deriving an 

exact expression for their partition function using Polychronakos's "freezing trick" . 
;> ■ From this expression we deduce several properties of the spectrum, and advance a 

''^ , number of conjectures that hold for a wide range of values of the spin M and the 

f^ I number of particles. In particular, we conjecture that the level density is Gaussian, 

^O ' and provide a heuristic derivation of general formulas for the mean and the standard 

^-p . deviation of the energy. 
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;h ' 1 Introduction 



The Haldane-Shastry (HS) chain [1,2] describes a fixed arrangement of equally 
spaced spin 1/2 particles in a circle with pairwise interactions inversely propor- 
tional to the square of the chord distance between the particles. The original 
interest of this model lies in the fact that the U —>■ oo limit of Gutzwiller's vari- 
ational wave function for the Hubbard model [3,4,5], which also coincides with 
the one- dimensional version of the resonating valence bond state introduced 
by Anderson [6], is an exact eigenfunction of the HS chain. The exact solv- 
ability of the HS chain was already proved in the original papers of Haldane 
and Shastry. A few years later Fowler and Minahan [7] used Polychronakos's 
exchange-operator formalism [8] to show that this model is also completely 
integrable. Although the obvious relation of the HS chain with the Sutherland 
(scalar) model of A^v type [9,10,11] was already remarked by Shastry, an ex- 
plicit quantitative connection was first established by Polychronakos through 
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the so-called "freezing trick" [12]. In the latter paper it is shown how to con- 
struct an integrable spin chain from a Calogero-Sutherland (CS) model of 
A]\f type with internal degrees of freedom ("spin") [13,14,15,16,17] by freez- 
ing the particles at the classical equilibrium positions of the scalar part of 
the CS potential. The first integrals of the spin chain are essentially obtained 
as the large coupling constant limit of the first integrals of the correspond- 
ing CS model. Polychronakos applied this technique to the original (rational) 
Calogero model of Ajsf type [18], constructing in this way a new integrable 
spin chain of HS type in which the spin sites were no longer equally spaced. 
In a subsequent publication [19], the same author gave a heuristic argument 
based on the freezing trick that relates the spectrum of the integrable spin 
chain with those of the corresponding scalar and spin dynamical models. 

Both the integrability and the spectrum of the Haldane-Shastry and Poly- 
chronakos spin chains can thus be obtained from the trigonometric and ra- 
tional CS spin dynamical models of An type. By contrast, the spin chains 
associated with the spin models of BCn type [20,21,22,23,24,25,26] have re- 
ceived comparatively little attention. This is in part due to the fact that, unlike 
their A^ counterparts, the BC^-type spin chains depend nontrivially on free 
parameters (one in the rational case and two in the trigonometric or hyperbolic 
cases). The integrability of the spin chain associated with the BCn rational CS 
model was established by Yamamoto and Tsuchiya [21] using the Dunkl oper- 
ator formalism [8,27,28], although, to the best of the authors' knowledge, the 
spectrum of this model has not been computed so far. The Haldane-Shastry 
(trigonometric) spin chain of BCn type was discussed by Bernard, Pasquier, 
and Serban [29], but only for spin 1/2 and with the assumption that the sites 
are equally spaced, which restricts the pair of free parameters in the model to 
just three particular values. Finkel et al. [25] recently discussed the integra- 
bility of the hyperbolic HS spin chain of BCn type, but did not examine its 
spectrum. 

In this paper we study the BCn version of the Haldane-Shastry spin chain for 
arbitrary values of the spin and the coupling constants. It turns out that there 
are actually four different BCn spin chains related to the original Haldane- 
Shastry chain, two of which are ferromagnetic and the other two antiferro- 
magnetic. We prove that these chains are exactly solvable provided that the 
sites are the coordinates of an equilibrium of a suitable scalar potential, which 
is the same for all four chains. In particular, for generic values of the coupling 
constants the sites are not equally spaced. In addition, we rigorously estab- 
lish the essential uniqueness of the equilibrium point of the scalar potential 
determining the chain sites. Using Polychronakos's freezing trick, we are able 
to derive an exact expression for the partition function of the models, thus 
establishing their solvability. From this expression, which is the main result of 
this paper, we deduce several interesting general properties of the spectrum. In 
the first place, the spectrum depends on the coupling constants only through 



their semisum (3, while for generic values of (3 the degeneracies of the energy 
levels depend only on the spin M and the number of particles N . Secondly, 
although the energy levels are in general unequally spaced, for f3 ^ N (and 
sufficiently large M) they cluster around an equally spaced set. In the third 
place, for half-integer spin the spectra of the two types of (anti)ferroniagnetic 
chains are exactly the same, even if their Hamiltonians differ by a nontrivial 
term. 

Apart from the rigorous results just mentioned, the evaluation of the parti- 
tion function for several values of M and A^ has led us to several conjectures 
regarding the spectrum. First of all, our calculations strongly suggest that the 
clustering of the levels around an equally spaced set when (3 ^ N occurs in 
fact for all values of the spin M. Secondly, even for moderately large values of 
N the level density follows a Gaussian distribution with great accuracy. This 
fact, which is the main conjecture of this paper, is reminiscent of the anal- 
ogous property of the "embedded Gaussian ensemble" (EGOE) in Random 
Matrix Theory [30]. It should be noted, however, that the essential require- 
ment defining the EGOE, namely that the ratio of the number of particles 
to the number of one-particle states tend to zero as both quantities tend to 
infinity does not hold in our case. If the level density is Gaussian to a very high 
degree of approximation (for sufficiently large N), it is fully characterized by 
the mean fi and standard deviation a of the energy. From a natural conjecture 
on the dependence of fi and a on the number of particles, we have derived 
general formulas expressing these parameters as functions of A^ and M. We 
have then checked that these formulas yield the exact values of ^ and a for 
a wide range of values of the spin and the number of particles. We have also 
rigorously proved (without making use of the previous conjectures) that the 
standard deviations for both types of (anti)ferromagnetic chains with integer 
spin exactly coincide, even if their spectra are essentially different. 

The paper is organized as follows. In Section 2 we introduce the Sutherland 
model of BCn type with internal degrees of freedom, and outline a proof 
of its integrability by expressing the Hamiltonian in terms of an appropriate 
commuting family of self-adjoint Dunkl operators. The spectrum of the latter 
model is determined in Section 3 by explicit triangularization of the Hamil- 
tonian. In particular, we compute the ground state energy in terms of the 
parameters of the model. The four types of Haldane-Shastry spin chains of 
BCn type, which are the main subject of this paper, are presented in Sec- 
tion 4. Section 5 is devoted to the calculation of the spectrum of the chains 
introduced in the previous section. We first provide a semi-rigorous detailed 
justification of the freezing trick, whose key points are the uniqueness of the 
equilibrium point of the associated scalar potential together with the knowl- 
edge of the full spectrum of the corresponding scalar and spin Sutherland 
models of BCn type. From the freezing trick we directly obtain an explicit 
expression for the ground state energy of the chains. We next make use of the 



freezing trick (which, by itself, does not completely determine the spectrum) 
to compute in closed form the partition functions of all four types of BC^ 
spin chains. In the last section we present concrete examples for spin 1/2 and 
1, which led us to formulate the general conjectures mentioned above. The 
paper ends with a technical appendix, in which we establish the uniqueness 
of the equilibrium point of the scalar potential determining the sites of the 
chains. 



2 The spin dynamical models 



In this section we shall study the integrability of the trigonometric Suther- 
land spin models of BC^ type. Each of these models describes a system of 
A^ identical particles with internal degrees of freedom ("spin") moving on a 
circle, subject to one- and two-body interactions depending on the particles' 
spatial and internal coordinates. We shall denote by S the finite-dimensional 
Hilbert space corresponding to the spin degrees of freedom spanned by the 
states |si, . . . , sn)i where —M < Si < M and M is a half-integer. We shall 
respectively denote by Sij and Si {i,j = 1, . . . , A^) the spin permutation and 
reversal operators, whose action on the basis of spin states is defined by 

i^ij Pl, . . . , Sj, . . . , Sj, . . . , Sjv/ pi, . . . , Sj, . . . , Sj, . . . , Sj\f) , , , 

Oj|Si, . . . , Sj, . . . , Si\f) = |Si, . . . , —Si, . . . , Sjv) . 

These operators are represented in S by (2M + 1) ^-dimensional Hermitian 
matrices. We shall denote by 6 the multiplicative group generated by the 
operators Sij and Si, which is isomorphic to the Weyl group of B^ type. We 
shall also use the customary notation Sij = SiSjSij. 

The -BCAT-type spin dynamical models we shall study in this section are col- 
lectively described by a Hamiltonian of the form 



H*^f = —^dl.+ a^ I sin ^ Xij (a — e Sij) + sin ^ xfj (a — e Sij 

i i^j 

+ b J2 sin'^Xi {b - e'Si) + b'J2 cos'^Xi ib' - t' S, 



(2) 
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where e, e' = ±1 are two independent signs, a, b, b' are real parameters greater 
than 1/2, and xfj = Xi±Xj. Here and in what follows, the sums (and products) 
run from 1 to A^ unless otherwise constrained. The potential in (2) possesses 
inverse-square type singularities at the hyperplanes Xi ± Xj = kn, Xi = k7i/2, 
with A; G Z. In fact, since the nature of these singularities makes it impossible 
for one particle to overtake another or to cross the singularities at Xi = kn /2, 
we can regard the particles as distinguishable and take as configuration space 



,N - ^ 



the set 

(7 = |x = {xi,...,xn) G 

The Hilbert space of the system may thus be taken as 7i = Ll{C) (X> S, where 



0<xi<---<xn<-\- (3) 



Ll{C) = \feL\C) 3 hm |x,±x,-A;7rr''|/(x)|, 3 hm |a;r'l/(x)| , 
3 hm \x, - 7r/2|-^'|/(x)| ; A; = 0, 1, 1 < ? ^ J < ivj. 

Note, in particular, that the physical wavefunctions vanish faster than the 
square root of the distance to the singular hyperplanes in their vicinity. 

Formally, the four Hamiltonians (2) can be represented as a single Hamiltonian 
H* = H*^, for an arbitrary choice of the signs e and e', provided that the 
parameters a, h and h' are also allowed to take negative values less than —1/2. 
We have preferred to use the more explicit representation (2) since, as we shall 
see in the following section, the spectrum of H^^, depends in an essential way 
on e and e'. It can be shown that the operator if*^, : Ti ^ Ti is equivalent to 
any of its extensions to spaces of symmetric or antisymmetric functions (with 
respect to both permutations and sign reversals) in L'l{C) ®S, where C is the 
A^-cube (— f , f )^ and Ll{C) is defined similarly to L^iC). We shall consider 
without loss of generality that H*^, acts in the Hilbert space 

7i,,. = K,.[Ll{C)®S), (4) 

where A^^' is the projection operator on states with parity e under simultaneous 
permutations of spatial coordinates and spins and e' under sign reversals. The 
latter operator is characterized by the relations 

KijK^^: = e S'jj A«/ , A'jA^e' = e'S'jA,,/ , (5) 

where Kij and Ki respectively denote the spatial coordinates' permutation 
and sign reversing operators, defined by 

\ ij J Jy'^ii • • • 1-^ii ■ ■ ■ -I -^j-i ■ ■ ■ 1 -^N j J y-^li . . . , 3Jj, . . . , Xj, . . . , Xjsjj , 

(i^i/)(xi,...,a;i,...,XAr) = /(a;i, . . . , -Xj, . . . ,a;Ar) . 



The relations (5) suggest the definition of a mapping *^; : 2) ® ^ — > 2) ® 6, 
where 2) denotes the algebra of scalar linear differential operators and .^ ~ 6 
is the multiplicative group generated by the operators Kij and Ki^ as follows: 

(dK,,,, ■ ■ ■ K,^,^Ki, ■■■KA\ = e'e"D Si^ ■ ■ ■ Si,S,^,^ ■ ■ ■ S,,,, , (6) 



where D E 1). This determines a hnear map A i-^ A* in 2) ® ^, which by 
Eq. (5) satisfies 

ak, = a:^,k,. (7) 

In particular, each of the physical Hamiltonians H*^, in (2) is the image under 
the corresponding star mapping of a single operator H , given by 

i? = — ^ (9^. + a ^ I sin"^ x^^ (a — Kij) + sin"^ x^- (a — Kij) 

' ' '^' , X (8) 

+ 6 ^ sin"^a;i {h- Ki) + h' Y^ cos^^x^ (6' - ii'i) . 



The integrability of the Hamiltonian (2) can be established by the same 
method applied in Ref. [25] to the hyperbolic version of //*_, based on the 
fact that H can be expressed as the sum of the squares of the commuting 
Dunkl operators 



^fc = i <9^fe + a ^ 



[l - i cot Xj^i) Km + (1 - i cot xti) Kki 

+ h {I - icot Xk) + h' {I + itdJiXk) Kk-la^Kki. (9) 



Kk 



These operators are related to the hyperbolic Dunkl operators J^ of Ref. [25] 
by Jfc(x) = — Jfc(ix). The commutativity of the Dunkl operators Jk implies 
that the operators 

h = llJk''^ p=l,...,7V, (10) 

k 

form a complete set of commuting integrals of motion of H = Ii. From this fact 
it follows, as in Ref. [25], that the corresponding operators {Ip)*^r, p = I, . . . ,N, 
act on the Hilbert space Tiee' and form a complete set of integrals of motion 
of the Hamiltonian H*^, = (h)*^,. 

To end this section, we shall prove that the integrals of motion {Ip)*^r are self- 
adjoint. Note first of all that, unlike the operators J^, the Dunkl operators (9) 
and hence the integrals of motion (10) are self-adjoint. Since, furthermore, Ip 
and Aee' are self-adjoint and commute with one another (in fact, Ip commutes 
with Kij and Ki by Lemma 4 of Ref. [25]), we have 



On the other hand, since {Ip)*^r also commutes with A^^r we obtain 

;(/,):,, A,,,) ^ = (a,,, (/,):,,) ^ = {IpyJa.,, 



from which it follows that (Ip)*^/ = {Ip)te' by Lemma 1 of Ref. [25]. 



3 Spectrum of the spin dynamical models 



In this section we sliall compute tlie spectrum of tlie trigonometric Suther- 
land spin models of BC^ type (2). The results of this section will be used in 
Section 5 to derive the asymptotic behavior of the partition function of these 
models in the large coupling constant limit. 

The computation of the spectrum of the Hamiltonian (2) is analogous to the 
corresponding computation for the hyperbolic model studied in Ref. [25], in 
spite of the fact that the boundary conditions are different. The starting point 
of this computation is the invariance under the Dunkl operators Jj of the 
finite-dimensional spaces 

T^fc = (0(x) exp (2i ^ n^Xj) Uj = -k,-k+l, . . . ,k , j = 1, . . . , iV^ , (11) 
j 

where 

0(x) = Y]^ I sinx^- sinx^l'' ■ Y[ \ sina;j|^| cosx^l^ , (12) 

i<j i 

for all nonnegative integer values of k. It follows that the operator H = Ii 
preserves the spaces 7^^ for all k. Since H commutes with A^^', Eq. (7) implies 
that 

H:^,[A,,,{ip\a))] = A,,,[iHip)\a)] , (13) 

for all ip G Lq{C) and \a) E S. Hence the Hamiltonian H*^, leaves invariant 
the infinite increasing sequence of finite-dimensional spaces 

Mk,ee' = Ke' ijlk ® S) , /c = 0, 1, . . . , (14) 

and is therefore exactly solvable in the sense of Turbiner [31,32]. 

We shall next construct a (non-orthonormal) basis B of the Hilbert space 
Lq(C) in which H is represented by a triangular infinite-dimensional matrix, 
thereby obtaining an exact formula for the spectrum of this operator. To this 
end, note that the (scaled) exponential monomials 

/n(x) = 0(x) exp (2i Y^ rijXj^ , n = (ni, . . . , un) , rij eZ, (15) 

j 

span a dense subspace of the Hilbert space Lq(C). We can introduce a partial 
ordering -< in the set of exponential monomials (15) as follows. Given a mul- 
tiindex n = (ni, . . . ,11^) G Z^, we define the nonnegative and nonincreasing 
multiindex [n] by 

[n] = (|nij,...,|ni^|), where In^J > ■ ■ ■ > |ni^| . (16) 

If n, n' G [Z^] are nonnegative and nonincreasing multiindices, we shall say 
that n -< n' if rii — n'l = ■ ■ ■ = nj_i — n',_i = and rii < n'. For two 



arbitrary multiindices n, n' G Z^, by definition n -< n' ii and only if [n] -< [n'\. 
Finally, we shall say that /„ -< /„/ if and only ii n -<n' . Note that the partial 
ordering -< is preserved by the action of the Weyl group ^, i.e., if /„ -< /„' 
then Wfn -< Wfn' for all W e ^. 

We can take as the basis B any ordering of the set of exponential monomi- 
als (15) compatible with the partial ordering -<. This follows from the fact 
that 

^/n = E%./n+ E ^nfn', n G Z^, (17) 

n' ^n 

where \[n],i and c" are real numbers (cf. Proposition 2 of Ref. [25]). The 
numbers \m,i ijn G [Z^]) are explicitly given by 

'2mi + h + h' + 2a(N + i + l- ij-{m.i) - 2i{mi)) , m^ > , 

Ami = < (18) 

-b-b' + 2a{i - N) , m^ = 



) 



where we have used the following notation: 

#(s) = card{i \mi = s} , i{s) = min{i \mi = s} . 

For instance, if m = (5, 2, 2, 1,1,1, 0) then #(1) = 3 and i{l) = 4. It will also 
be convenient in what follows to take i{s) = +oo if m, 7^ s for alH = 1, . . . , A^. 
Equation (17) implies that the operator H is represented in the basis B by an 
upper triangular matrix with diagonal elements 



-E„ = E-^ 



2 

n],i 



From the previous formula it is straightforward to deduce the following more 
compact expression for the eigenvalues En of the operator H: 



E„ = E 2[n], + 6 + 6' + 2a(iV-i) . (19) 



2 



Indeed, ii m = [n] G [Z^] and rrik-i > rrik = ■ ■ ■ = irik^p > rrik^p^i > then 
i{mk+j) = k and ^{rrik+j) = p + 1 for j = 0, . . . ,p, so that 

Kn,k+j = 2mk+j+b+b'+2a{N-k-p+j) = 2mk+p-j+b+b'+2a(^N-{k+p-j)'^ 
and hence 

k+p k+p 

J2\l^^ = J2{2m, + b + b' + 2a{N-t)) . (20) 

i=k i=k 

If, on the other hand, ruk-i > ruk = ■ ■ ■ = m^ = 0, Eq. (20) follows directly 
from (18). This completes the proof of the formula (19). 



Let us now compute the spectrum of the Hamiltonian (2) in Tiet'- Note, first 
of all, that the states of the form 

A«'(/n|si,...,s^)), ne[Z^], (21) 

span a dense subset of the Hilbert space Ti^e', by the analogous property of the 
functions (15) in Lg(C). The states (21), however, are not linearly independent 
(in particular, some of them vanish if e or e' are negative). A (non-orthonormal) 
basis of Tiee' may be obtained from the states (21) by imposing the following 
conditions on the spin vector |si, . . . , sn) (cf. Proposition 3 of Ref. [25]): 

i) Si — Sj > 6-i^e , if ni = rij and i < j ', (22a) 

ii) s, > i 5_i,,, , if n^ = , (22b) 

where 6 is Kronecker's delta. Indeed, if n^ = rij with i < j we can clearly 
permute the i-th and j-th particles (if necessary) so that Sj > Sj, leaving the 
state (21) invariant up to a sign. If, in addition, e = — 1 we must have Sj > Sj 
by antisymmetry under permutations. Likewise, if n^ = we can assume that 
Sj > after a possible reversal of the sign of the coordinates of the i-th particle, 
which again preserves the state (21) up to a sign. Moreover, e' = —1 forces 
Sj > by antisymmetry under sign reversals. Note that when e = — 1, i.e., 
when the basis states (21) are antisymmetric with respect to permutations, 
the first condition implies the following restriction on the multiindex n G [Z^] : 

f2M + l, if ni>0 

*ini)<{ , ' ., ^ (23) 

[ M,,, if n, = 0, 

where M+ = [MJ + 1 and M_ = [M] . Here \_x\ and [x] denote respectively 
the integer part of x and the smallest integer greater than or equal to x. Let 
Bee' be any ordering of the set of states (21)-(23) compatible with the partial 
ordering -<. It follows from Eqs. (13) and (17) that the matrix of H*^, with 
respect to the basis Bee' is upper triangular, with eigenvalues given by 

E*^,{n;s) = J2{'^n, + b + b' + 2a{N-i))\ n e [Z^] , s = (si, . . . , sjv) • 

(24) 
It is worth mentioning at this point that, although a cursory inspection of 
the previous equation may suggest that the models (2) are isospectral, this is 
in general not the case. Indeed, condition (23) implies that many eigenvalues 
of the models with e = 1 are absent from the spectrum of the models with 
e = — 1. Besides, for a fixed value of e, by condition (22b) the degeneracy of 
the eigenvalues also depends on e' when M is an integer (see Eqs. (50) and 
(51) in Section 5 for the minimum degeneracy of each level). 

Since E*^,[n; s) is an increasing function of the components of the multiindex 
n, the ground state of the system is obtained when each component n^ takes 



the lowest possible value. Thus for e = l, ore = — 1 and A^ < M^/, we have 
n = 0. On the other hand, when e = — 1 and A^ > M^r condition (23) implies 
that 



2M+1 



2M+1 



M,, 



n 



mo, . . . , mo, mo - 1, . . . , mo - 1, . . . , 1, . . . , 1, 0, . . . , 



(25) 



where N = M,> + (mo - 1)(2M + 1) + r with r = 1, . . . , 2M + 1. The ground 
state energy E*^,^^^ is easily computed in this case using Eq. (24) with the 
multiindex n given in (25). We thus obtain 



K^min = I «'^' - 2aciV2 + i (3c2 - a')N + ^ /.mo 



4mQ(l — an) 



+ 6cmo + a/t + 2 + 2mop c + mo(l — clk,) — t:0'P 



(26) 



where c = a — b — b' — 2mo and 

K = 2M + 1 , 



e', M = 0,1,... 

2'2' 



(27) 



It can be easily shown that Eq. (26) is also valid when n = if we take mo = 0. 
Thus Eq. (26) yields the ground state energy in all cases provided that mo is 
defined by 

^N -My 



mo = 5, 



e,-i 



2M+1 



28 



4 The spin chains 



The Hamiltonian of the HS spin chains of BC^ type associated with the spin 
dynamical models (2) discussed in the previous sections is given by 

+ J2{(3 sin-^ e, + P' cos-' e.) (1 - eS,) , (29) 



where (3 and (3' are positive real parameters, ^^ = ^i±Q, and ^ = {^i, ■ ■ ■ , ^n) 
is the unique equilibrium point in the set C of the classical potential 

f/(x) = ^(sin-^ x^j + sin"2 x,^) + J^iP^ sin"' Xi + /?" cos"' Xi) . (30) 

i^j i 

It is important to note that the classical potential (30) is independent of e and 
e', and therefore the sites of the four chains (29) are the same. The existence of 



10 



a minimum of f/ in C for all values of (3 and (3' is a consequence of the positivity 
and continuity of f/ in C and the fact that it tends to infinity at the boundary 
of this set. The uniqueness of this minimum is proved in Appendix A. Note 
that, in contrast, the corresponding potential for the hyperbolic spin chain of 
BCn type treated in [25] admits an equilibrium point only for a certain range 
of values of 13 and (3' . 

The chains (29) with e = — 1 (respectively e = 1) are of antiferromagnetic 
(respectively ferromagnetic) type. Note also that the spin chain Hamiltonians 
^-e-e' and —\\ee' are related by 

h_,,_,, = -h,,, + 2\/(^), (31) 

where 

y (x) = ^ (^sin~^ x^^ + sin"2 xf^ + /? H sin^^Xj + /3' ^ cos^^x^ . (32) 

i^j i i 

On the other hand, for a fixed e the two chains h^ j- are essentially differ- 
ent. From Eq. (29) it immediately follows that the eigenvalues e^e/j- of the 
Hamiltonian h^^/ are nonnegative. Moreover, in the ferromagnetic case (e = 1) 
the ground state energy clearly vanishes, since states symmetric under per- 
mutations with parity e' under spin reversals are annihilated by the Hamilto- 
nian. Equation (31) implies that the maximum energy of the antiferromagnetic 
chains is 2V{$,). 

The integrability of the spin chain (29) with M = 1/2 and e = e' = — 1 
was proved in Ref. [29] only for the special values (3/2, 1/2), (3/2,3/2), and 
(1/2, 1/2) of the pair {P,P'), for which the corresponding sites 

in in / 1\ tt 

«-- -T-; 1 = 1 iV, 



2iV + l ' 2A^ + 2 ' V 2/ 2A^ 



) -^ ' ) 



are equally spaced, as in the original Haldane-Shastry chain [1,2]. As we shall 
see below, the discussion of the integrability of the HS spin chains of BCn 
type (29) for arbitrary values of the parameters f3 and /?' is completely analo- 
gous to that of Ref. [25] for the hyperbolic version of h 

Let us begin by defining the operators Jj {i = 1, . . . ,N) and Ip {p G N) by 

Ji = i <9^, + aJi , \p = Yl ■^?- 

i 

The operators (Ip)*^/ clearly commute with one another, since [(Ip)*^/, (Ig)*^/] 
is the coefficient of a^^^'^'^^ in the expansion in powers of a of the identity 
[(/p)*^,, (/g)*e'] = 0. Hence the operators 

(w:° = (ip):.'U > peN, (33) 



11 



also commute with one another. We shall now prove that the operators (33) 
form a commuting family of integrals of motion for the spin chain Hamiltonian 
hge'. Note that this result does not follow trivially from the previous assertions 
since, in contrast with the dynamical case, (li)*^/ = U{^) is a constant and 
therefore does not coincide with he^'. 

The starting point in the proof of the commutativity of h^e' and (lp)*°/ is the 
following expansion of H in powers of a: 



H = -j:dl-aH + a'U{^) 



(34) 



where 



H = E 



sin ^ Xj^j Kij + sin ^ xfj Kij 



+ ^(/5sin-2a;i + /5'cos"2^i)fs:i, (35) 



U is defined in Eq. (30), and we have set 



a a 



(36) 



Note that h^e' = — H*° + V{$,), so that the integrability of h^e' follows from 
that of H*^,. Arguing as in Ref. [25] it is straightforward to show that 

dU 

i ^ 

for certain operators Cp^j in 2D (g) .^. Setting x = ^ in the previous identity, 
it follows that (Ip)*^/ commutes with h^^/. Note finally that the first integrals 
(lp)*° are clearly self-adjoint, since they are equal to the coefficient of a^^ in 
the corresponding self-adjoint operators (/p)*^/ evaluated at the equilibrium 
point $,. 



5 Partition function and spectrum of the spin chains 



In this section we shall compute the partition function of the HS spin chains of 
BCn type (29) by using Polychronakos's freezing trick [12,19] applied to the 
spin dynamical models discussed in the previous sections. We shall first provide 
a detailed heuristic justification of the freezing trick in the present context. 
Our calculation relies on the computation of the large coupling constant limit 
of the partition functions of the spin dynamical models (2) and the scalar 
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Sutherland model of BCn type 



sm ^ Xjj- + sin ^ x^^ 



+ h{h -l)Y. sin'^x, + h'{h' -l)Y. cos-^x, 



acting on the Hilbert space Ll{C). Using the definition (36) of /3 and 13' we 
obtain 

^s = -E^',+«'f^W-«^W' (38) 

i 

with f/(x) and V^(x) respectively given by Eqs. (30) and (32). From Eqs. (34) 
and (38) it follows that 

K, = -J2dl-a H:, + a^f/(x) = H, + a(\/(x) - H:,,) , (39) 

i 

where H*^, is assumed to act in the Hilbert space 7i = Lq{C) C?> S. Let 
|-?/'j(x) \ be a basis of eigenfunctions of H^, and let I |o'm',j) [ . , with d = 
(2M + 1) , be a basis of eigenfunctions of h^e', so that 

The set <^-?/'j(x)|o"„' ,) > is thus a basis of the Hilbert space Ti., and 

ifs(V^i(x)|o-,e'j)) = Ei^i{x.)\a,,>j) , 

since i/g does not act on the spin variables. 

From Appendix A it follows that the classical potential a^ t^(x) — a V(x) has 
a unique equilibrium point (actually, a minimum) x(a) in the set C provided 
that a > max(l//9, 1//3', 1). The freezing trick is based on the fact that for 
a ^ 1 the eigenfunctions of Hs are all sharply peaked around the equilibrium 
x(a)- Since x(a) = ^ + 0(a~^], for a ^ 1 we have 



\/(x) - H:,, U,(x)|a,,,,)) = ^,(x) \/(x) - H:,, |a,,,,) 



~ V^i(x)[\/(^) - H*°,]|(T,,/j) = V^i(x)(h,,.|a,e'j)) = e,e'j^i(x)|o-,,/j) . 

Thus for a ^ 1 the Hamiltonian H*^, is approximately diagonal in the basis 
|-?/'j(x)|o"ee',j) [. . , with eigenvalues approximately given by 

E*^,^i^ ~ Ei + ae«'j , ieN, l<j<d, a>l. (40) 

Taking into account that Ce^'j is independent of a, we immediately obtain the 
following exact expression for the eigenvalues e^^'j of the spin chain Hamilto- 
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nian (29): 

e,,,, = lim - {El,,^ - Ei) . (41) 

Using Eq. (40), we can easily derive the ground state energy eee',min of the spin 
chains (29). Indeed, clearly the ground state energy -E*e'min given by Eq. (26) 
is achieved when both Ei and e^e',j in Eq. (40) attain their minimum values 
-Emin and eee'Mn- From Eqs. (24) and (36) it follows that both -E*e',jj and E^ 
are polynomials of the second degree in a with the same leading coefficient. 

Hence 

1 

ee.',min = jim - (^;,,„,i„ - ^min) • (42) 

Since -Emin is obtained by setting n = in Eq. (24), the coefficient of a in -Emin 
vanishes. From the previous equation, it follows that eee',nnn is the coefficient 
of a in -E*^/jnm5 namely (cf. Eq. (26)) 



^ee'min ''^0 



AN^ + 4(2/3 -l)N +'^(k- 2mo(6/5 + 2moK - 3) 



-2p(2/3 + mofi:-l)-p'J, (43) 
where k, p and mo are defined in Eqs. (27) and (28), and we have set 

'^=\{P + P'). (44) 

Note that (as remarked in the previous section) the ferromagnetic ground state 
energy vanishes, since ttiq = when e = 1. 

We emphasize that Eq. (41) cannot be used directly to compute in full the 
spectrum of h^e', since it is not clear a priori which eigenvalues of H*^, and Hg 
can be combined to yield an eigenvalue of h^^/. The importance of Eq. (41) lies 
on the fact that it can be used as the starting point for the exact computa- 
tion of the partition function of the spin chain h^e/, which in turn completely 
determines the spectrum. 

Let us denote by Z^, Z*^, and Z^^/ the partition functions of the scalar Suther- 
land Hamiltonian (37), the Sutherland spin dynamical model (2), and the spin 
chain Hamiltonian (29), respectively. From Eq. (40) it follows that Z*^,[T) ~ 
Zs(T) Zee'(T/a), and hence 

Z.AT) = jun^f^ . (45) 



a— >oo 



Recall that the spectrum of the Hamiltonian H*^, of the spin dynamical model 
is given by Eq. (24), where n G [Z^] is a nonnegative nonincreasing multiindex 
(satisfying conditions (23) if e = —1), and s = (si, . . . , stv), —M < Sj < 
M, satisfies (22a)-(22b). The leading terms in the expansion of E*^,[n; s) are 
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therefore 

E*^,{n; s) - a^Eo + SaJ^niiP + N - i) , (46) 

i 

where Eq = 4:J2i{f3 + N — i^ is a constant independent of n. The eigenvalues 
E{n) of the scalar Sutherland Hamiltonian Hg are given by the right-hand side 
of Eq. (24), where the multiindex n G [Z^] is now unrestricted. Thus E{n) 
also satisfies Eq. (46) for a ^ 1. 

Let us start by computing the large a limit of the denominator in Eq. (45). 

Note, first of all, that the constant Eq can be dropped from both E*^, (n; s) and 

E{n) without affecting the value of Z^^i{T). Using the asymptotic expansion 

of E{n) and setting 

q = e-8/(^BT) (47) 

we immediately obtain 

Z,{aT):^ Y, q^ . (48) 

Defining Pi = rii — Uj+i, 1 < i < N — 1, and pjy = ''^Af we have 

« i<j 3 3 

and hence 



Pi,...,Pjv>0 J i Pi>Q 



(49) 



Let us compute next the partition function Z*^,{aT) of the Sutherland spin 
dynamical model (2) for a ^ 1. To this end, it is convenient to represent the 
multiindex n G [Z^] appearing in Eq. (46) as 



fcl fc2 



/, " ^ / " ^ / " ^\ 

n = (^mi, . . . ,mi,m2, . . . ,m2, . . . ,mr, . . . ,mr j, (50) 

where mi > 1712 > ■ ■ ■ > rrir > 0, and ki = i^{mi) G N satisfies ki + - ■ ■+kr = N 
(together with condition (23), if e = — 1). Thus 

E:,,(n;s)~8a5]m, J2 (P + N - j) 

*" /_ Ik- *~^ \ *" 

= 8a E rriiki {P + N---^~YkA =80^ rriiPi . 
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Let k = {ki, . . . ,kr), and denote by dee'{k, rrir) the cardinal of the set of spin 
quantum numbers s satisfying conditions (22a)-(22b) for the muhiindex (50), 
namely 



d,,(k,mr.) = (^M+I+JmC^i-D) . . . (2M+l+5,„(fc.-l) 



nir > 0] 



(51a) 



d,4k, 0) = (^2^'^^+i+^^M('^i-i)j . . . [''''^'+l';['^-'-'>) ['''^'^X'^-'^) ■ (51b) 



\ ki J ' ' ' \ kr^l J \ 

The partition function Z*^,{aT) is therefore given by 



keVN mi>-->mr>0 



i=l 



r-1 



E E d,Ak,mr.)l[r-''+ i: E d,Ak,0)l[q^^''% 

fegPjv m,i>--->mr>0 i=l fegPjv m,i>--->mr-i>0 



j=l 



where we have denoted by Vn the set of partitions of the positive integer A^. 
Since 

s 

E nr'"= E n? '- = E nn/^" 

mi>-->ms>0 j=l pi,...,ps>Oi=l pi,...,ps>0 i=l j=i 



P^E'^' ^ g^^- 



pi,...,Ps>Ojr = l j = lp^>0 j = l 



g^. ' 



where 






(52) 



using Eqs. (51) we finally obtain 



z:AaT) ^ E 

{kl,...,kr)€PN 



M^/+5l,E(fer- 



1)\ , (2M + l + 5l,,(kr-l)\ V 

y ^ V kr ) l_ (jNr 



r-1 



2M+l+5i,,(fcj-l)\ Q^ 



Af, 



1-g 



N, 



■ (53) 



Equations (45), (49) and (53) yield the following exact formula for the partition 
function of the HS spin chain (29): 



N 



z«'(T) = n 



i=l 



+ 



l-q 



i\l3+N-\{i+l) 



2M+l+5i,^(kr-l)\ ^ 



l-qNr 



E 

J (ki,...,kr)&VN 

r-1 

n 

J j=l L 



M^,+Si^c{kr-l) 



N, 



2M+l+5i,,{fcj-l)\ ([ 



N, 



1-q 



N, 



■ (54) 



From the previous formula, which is in fact the main result of this paper, one 
can infer several remarkable properties of the spectrum of the spin chain (29) 
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that we shall now discuss. First of all, for half-integer M the partition func- 
tion (54) does not depend on e', since in this case M± = M + 1/2. Hence the 
spectrum of the spin chain (29) is independent of e' when M is a half-integer, 
a property that is not immediately apparent from the expression of the Hamil- 
tonian (29). Secondly, all the denominators 1 — g^*, 1 < /c < r, appearing in 
the second line of Eq. (54) are included as factors in the product in the first 
line. Hence the partition function (54) can be rewritten as 

ZAT)= E 4e',^(M)g^^ (55) 

where Ss is given by 

N _ 1 

e5 = Y.^^^{(^ + N--{t + l)), (5 = (<5i, . . . , (5^) , (56) 

and the degeneracy factor ^^^/^(M) is a polynomial of degree N in M. There- 
fore, for all values of e, e' and M , the spectrum of the spin chain is contained 
in the set of 2^ numbers Se^, 5 G {0, 1}^. Moreover, for generic (sufficiently 
large) values of the spin M, the spectrum exactly coincides with the above set 
of numbers, the values of e, e' and M affecting only the degeneracy d^^i^s{.M) of 
each level. In particular, from the previous observation and Eqs. (31) and (56) 
we immediately obtain the following exact expression for the constant V{$,) 
(i.e., half the maximum energy of the antiferromagnetic chain): 

V{i) = AY^i{p + N - -{i + l)) = -N{N +1){2N + ?,p -2) . (57) 

From Eq. (56) it follows that the energies of the spin chains (29) are of the 
form 8(j/9 + A;), with j, /c nonnegative integers. Since the coefficients of the 
powers of q appearing in Eq. (53) are independent of /5, it follows that for 
generic^ values of (3 the degeneracy of the levels depends only on the spin M 
and the number of particles N . 



6 Discussion and conjectures 



In this section we shall present several concrete examples in which we apply the 
formula (54) to compute the spectrum of the spin chains (29) for certain values 
of A^ and M . These examples strongly suggest a number of conjectures that 
shall be discussed in detail at the end of this section. We shall restrict ourselves 



^ More precisely, for all real values of /3 except for a finite (possibly empty) set of 
rationals. 
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to the antiferroinagnetic chains h_^±, the properties of their ferromagnetic 
counterparts following easily from the relation (31). 

Example 1 . The structure of Eq. (54) makes it straightforward to compute 
the spectrum of the spin chains for any fixed number of particles as a function 
of the spin. For instance, for A^ = 3 sites and integer M the energies (divided by 
8) of the spin chain h__ are 0, ]3+2, 2]3+?>, 'i]3+?>, 3^5+5, A]3+h, 5^5+6, G^^+S, 
with respective degeneracies 

iM(M-l)(M-2), |M(M2-1), |M(M + 1)(11M-2), 

\ M{M + 1)(7M - 4), i M{M + 1)(7M + 11), i M{M + 1)(11M + 13), 

I M{M + 1)(M + 2), |(M + 1)(M + 2)(M + 3). 

Note that in this case all energy levels Sss, with ss given by (56), are attained 
for M > 3, in agreement with the general discussion of the previous section. 
For M = 1 and a few values of [3 and j3', we have numerically computed the 
spectrum of the spin chain (29) by representing the operators Sij and Si as 
27 X 27 matrices. The results obtained are in complete agreement with those 
hsted above. 

For a fixed value of the spin M, we have not been able to find an explicit for- 
mula expressing the energies and their degeneracies as functions of the number 
of particles A^. However, if we fix M the spectrum can be straightforwardly 
computed from Eq. (54) for any given value of N. We shall next present two 
concrete examples for the cases M = 1/2 and M = 1. 

Example 2: spin 1/2. In this case we have computed the partition func- 
tion Z_ -I- for up to 20 particles (recall that for half-integer M the chains with 
e' = ±1 have the same spectrum). For instance, for N = 6 the antiferromag- 
netic spin chain energies (divided by 8) and their corresponding degeneracies 
(denoted by subindices) are given by 

(9;5 + 32)2, (10;5 + 36)2, (li;5 + 38)2, (li;^+41)4, (12;5 + 38)i, (12;5 + 43)6, 
(13^+43)3, (13:^ + 46)6, (14:^ + 46)4, {14^ + 50)^, (15;5 + 47)2, {15p + 50^, 
(15^+55)6, (16;5 + 51)2, (16^+55)5, (17^ + 53)i, (17^+56)4, (18^ + 58)3, 
(19^+61)2, (20;5 + 65)i, (2i;5+70)i. 

The number of levels increases rapidly with the number of particles A^. For 
example, if A^ = 10 the number of levels (for generic values of /3) is 136, 
while for A^ = 20 this number becomes 7756. It is therefore convenient to 
plot the energy levels e^ and their corresponding degeneracies di, as is done 
in Fig. 1 for A^ = 10 particles. Note that Eq. (56) implies that when (3 ^ N 
the levels cluster around integer multiples of 8/9. In fact, for all A^ up to 
20 we have observed that these integers take all values in a certain range 
io) Jo + 1) • • • ) A^(A^ + l)/2; for example, in the case A^ = 6 presented above 



Jo = 9. 




di 



Fig. 1. Energy levels e^ and degeneracies di of the antiferromagnetic spin 1/2 
chain h_^-|- for iV = 10 particles and j3 = v2- 



Example 3: spin 1. We have computed the partition functions 7._ ± of the 
spin chains h_ j- with spin M = 1 for up to 15 particles. As remarked in the 
previous section, for integer M the partition functions Z_ ± are expected to be 
essentially different. This is immediately apparent from Fig. 2, where we have 
compared graphically the energy spectra of the even and odd spin chains Z_ -t 
with (3 = v2 for A^ = 10 particles. However, we shall prove in what follows 
that the standard deviation of the energy is exactly the same for both chains. 
This rather unexpected result will be relevant in the ensuing discussion of 
the level density (see Conjecture 2 below). We also note that, just as for spin 
1/2, for A^ up to (at least) 15 and /5 ^ A^ the energy levels cluster around an 
equally spaced set of nonnegative integer multiples of 8/3. 

The previous examples for spin 1/2 and 1 suggest several conjectures on the 
spectrum of the (antiferromagnetic) HS spin chains of BCn type that we shall 
now present and discuss in detail. 

Conjecture 1 . For j3 ^ N, the energies cluster around an equally spaced 
set of levels of the form 8j^, with j = jo, jo + !,•••, A^(A^ + l)/2. 

In fact, for sufficiently large values of the spin M this assertion (with jo = 0) 
follows directly from Eq. (56). Our calculations for a wide range of values of 
A^ and M fully corroborate the above conjecture. 

Conjecture 2. For N ^ 1, the level density follows a Gaussian distribution. 

More precisely, the number of levels (counting their degeneracies) in an interval 
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di 



Fig. 2. Comparison of the energy levels e, and degeneracies di of the antifer- 
romagnetic spin 1 chains h (left) and h |- (right) for A^ = 10 particles and 

;5 = \/2. 



I is approximately given by 



{2M+lf jAf{e-fi,a)de, 



(5J 



where 



^{e;i2, 



a 



1 (e-M)^ 

e ~^^^ 



aV27i 

is the normal (Gaussian) distribution with parameters /x and a respectively 
equal to the mean and standard deviation of the energy spectrum of the spin 
chain. Although the shape of the plots in Figs. 1 and 2 make this conjec- 
ture quite plausible, for its precise verification it is preferable to compare the 
distribution function 



(60) 



F^(e) = / Afit; fi, a) dt 

of the Gaussian probability density with its discrete analogue 

F(e) = (2M+l)-^ Y: di, 

i; ei<e 



(61) 



where di denotes the degeneracy of the energy level Cj. Indeed, our computa- 
tions for a wide range of values of M and A^ > 10 are in total agreement with 
the latter conjecture for all four chains (29). This is apparent, for instance, in 
the case /3 = \/2, M = 1/2, and iV = 10 presented in Fig. 3. The agreement 
between the distribution functions (60) and (61) improves dramatically as A^ 
increases. In fact, their plots are virtually undistinguishable for A^ > 15. 
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F^ie), F(e) 




1750 2000 2250 2500 2750 3000 3250 



Fig. 3. Distribution functions F_\f{e) (continuous line) and F{e) (at its disconti- 
nuity points) for ]3 = V2, M = 1/2, and N = 10. 



It is well known in this respect that a Gaussian level density is a characteristic 
feature of the "embedded Gaussian ensemble" (EGOE) in Random Matrix 
Theory [30]. It should be noted, however, that the EGOE applies to a system 
of N particles with up to n-body interactions (n < N) in the high dilution 
regime N -^ oo, k, ^ oo and N/n -^ 0, where n is the number of one-particle 
states. Since in our case k, = 2M + 1 is fixed, the fact that the level density is 
Gaussian does not follow from the above general result. A study of the energy 
spectrum of the spin chains (29) in the framework of Random Matrix Theory 
is nonetheless worth undertaking, and will be the subject of a subsequent 
publication. 

If Conjecture 2 is true, the level density for large N is completely characterized 
by the parameters fi and a through the Gaussian law (59). It is therefore of 
great interest to compute these parameters in closed form as functions of A^ 
and M. To this end, let us write 



,± 



E 



hail + Sii) +hiJl + Si. 



Hj 



^Ij 



+ Y.h,{lTSi, 



(62) 



where the constants hij, hij and hi can be easily read off from Eq. (29). We shall 

begin by computing the average energy fi = /i_ of the odd antiferromagnetic 

spin chain h for integer spin. Using the formulas for the traces of the spin 
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operators given in Table 1 we immediately obtain 

2(M + 1 



/i_ = (2M + l)-^trh__ 
2(M + 1 



2M+1 



Vi^) 



2M + 1 
4(M + 1 






3(2M+1) 



N{N+l){2N + ?>(3-2), 



M en, 

(63) 



where we have used the explicit expression (57) for V{$,). On the other hand, 

Table 1 

Traces of products of the spin operators. 



Operator 




Trace (integer M) 




Trace (half-integer M) 


S^ 




(2M + 1)^-1 







Jij ) Jij 




(2M + 1)^-1 




(2M + 1)^-1 


biOj 




(2M + 1)^-2+25., 




(2M + l)^<5i,- 


^ij^kj ^ij 


Sk 


(2M + 1)^-2 







Si j Ski 




(2M + 1)^-2 




{2M + l)''-^{l-5^k5,i){l-6u6jk) 


SijSkl, Si 


jSki 


(2M + 1)^^-2+25^^5,^ 


+2<5ii<5jfc 


{2M + l)N-2+25,k&ji+2Sa&jk 



the average energy /x__|_ = yU+ of the even chain h_+ is given by 

;.+ = (2M + l)-^trh_+ = ?i^^E(^., + M + ^^E^^ 



2M + 



«7^i 



2M+1 



(M+l)\/(^)-Si , MgN, 



(64) 



where Si = X]j ^i is obviously independent of the spin M . Similarly, for half- 
integer spin the formulas for the traces of the spin operators in Table 1 yield 
the following expression for the mean energy /i_,± = /x±: 



/i± 



2M+ 1 



2(M+l)\/(^)-Si 



2' 2' 



(65) 



Let us turn now to the (squared) standard deviation of the energy, given by 

tr(h2 J 



a_ I = cr. 



(trh„,±)2 



'±~^ (2M+1)^ (2M + 1)2^- 
For integer spin, a long but straightforward calculation using the formulas in 
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Table 1 yields 

4M(M + 1) 



^± 



{2M + iy 



i^j i 



4M(M + 1)^ _ _ 



(2M + 1) 

(66) 
Since S2 does not depend on M, the above equation completely determines 
the dependence of a± on the spin. An important consequence of the previous 
formula is the equality of the standard deviation of the energy for the even and 
odd antiferromagnetic chains (for half-integer spin, this follows trivially from 
the fact that the even and odd chains have the same spectrum). This result is 
quite surprising, since for integer spin the energy spectra of the chains h_^-|- are 
essentially different, cf. Fig. 2. For half-integer spin, an analogous calculation 
yields the expression 



cr_i 



4M(M + r 
(2M + 1)2 



M(M + 1) 



M=i.5..... (67) 



where 



7 E ^*^ - E ^ijh (6^ 



4 



is independent of the spin. As before, Eq. (67) fixes the dependence of (T-j- on 
the spin. 

We still need to evaluate Si(A^), S2(N) and H^iN) in order to determine the 
dependence on A^ of /i± and a± in all cases. Although we have not been able 
to compute these quantities in closed form, in view of Eq. (63) it is natural to 
formulate the following conjecture: 

Conjecture 3. The average energy n± and its squared standard deviation 
(T^ depend polynomially on N. 

In fact, since e_^-|-;max = 2V^(^) is a polynomial of degree 3 in A^ by Eq. (57), it 
follows that the degrees in A^ of fi± and o"^ cannot exceed 3 and 6, respectively. 
The latter conjecture and this fact allow us to determine the quantities Sj(A^) 
by evaluating fi± and o"^ for A^ = 2, . . . , 8 and M = 1/2, 1 using the exact 
formula (54) for the partition function (cf. Eqs. (64), (66) and (67)). The final 
result is 

El = E3 = 2N{2^ + N -1), 



S 



2 



4A^ 
"9~ 



>^2 



2(2A^' + 3N + 13)/? + (A^ - 1)(5A^' + 7N + 20) f3 (69) 

+ -(N - 1)(8N^ + 3N^ + 13N - 12) 
5 

These expressions, together with Eqs. (63)-(67), completely determine fi± and 
a± for all values of M and A^. We have checked that the resulting formulas 
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yield the exact values of fi± and (x± computed from the partition function (54) 
for a wide range of values of M and A^. This provides a very solid confirma- 
tion of Conjecture 3. Let us mention, in closing, that formulas analogous to 
(63)-(67) expressing the mean and standard deviations of the energy for the 
ferromagnetic chains h+^± can be immediately deduced from the previous ex- 
pressions and Eq. (31). 
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A Uniqueness of the equilibrium of the classical potential 



In this Appendix we shall prove that the classical potential (30) has exactly 
one equilibrium point in the set C. We have already seen in Section 4 that U 
has at least one minimum in C. We shall now prove that the Hessian of U is 
positive-definite in C, which implies that all the critical points of f/ in C must 
be minima. This implies that U has exactly one critical point (a minimum) in 
the set C. 

If f(t) = sin~^ t, we can express the second partial derivatives of U as follows 
'^' (A.l) 






/"(4) - n^r 



Note that f"{t) = 2 csc^ t(l + 3 cot^ t) is strictly positive for all values of t, and 
therefore {d'^U) / {dxj) > for all i. By Gerschgorin's theorem [33, 15.814], the 
eigenvalues of the Hessian of U lie in the union of the intervals 



dx'i 



-li 



dx'i 



+ li 



, where 7i = E 






1 = 1,. ..,N. 



Since 



^ -%> p'f'ix,) + p"f"{^ - x;) > 0, 
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all the eigenvalues of the Hessian of U are strictly positive. This establishes 
our claim. 
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